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We analyse the prospect of extending the reach for squarks and ganginos via associated production 
at a ^/s = 100 TeV proton-proton collider, given 3 ab“^ integrated luminosity. Depending on the 
gluino mass, the discovery reach for squarks in associated production with a gluino can be up to 37 
TeV for compressed spectra (small gluino-LSP mass splitting), and up to 32 TeV for non-compressed 
spectra. The discovery reach for Winos can be up to between 3.5 and 6 TeV depending on squark 
masses and Wino decay kinematics. Binos of up to 1.7 TeV could similarly be discovered. Squark- 
gaugino associated production could prove to be the discovery mode for supersymmetry at a 100 
TeV collider in a large region of parameter space. 


I. INTRODUCTION 

Observational evidence for low energy Supersym¬ 
metry (SUSY) remains elusive. Current LHC data 
constrains strongly-interacting superpartner masses 
to lie near or above a TeV, disfavoring electroweak- 
scale SUSY in a wide variety of models. It is there¬ 
fore becoming increasingly well motivated to con¬ 
sider the possibility that the superpartner masses lie 
above ~ 1 TeV, perhaps evading the kinematic reach 
of LHC-14. This has prompted numerous studies of 
the SUSY discovery potential of future hadron col¬ 
liders, which have demonstrated that a ^/s = 100 
TeV collider can extend the kinematic reach for su¬ 
perpartners into the multi-TeV range [1-11]. 

Previous studies of SUSY at future hadron collid¬ 
ers have focused primarily on pair production, either 
of colored superpartners [1, 2, 4] or of electroweak- 
inos [3, 6-10]. In this paper, we instead examine the 
reach of a y^s = 100 TeV collider for associated pro¬ 
duction of a heavy squark along with a lighter gaug- 
ino. This production channel is particularly note¬ 
worthy if the squark masses are (!l(10)’s of TeV, such 
that squark pair production is kinematically inacces¬ 
sible at i/s = 100 TeV. Spectra where squarks are 
hierarchically heavier than the gluino/electroweak- 
inos are predicted in many SUSY breaking models 
such as anomaly mediation [12, 13] or more gen¬ 
eral “mini-split”-type scenarios [14-16]. Moreover, 
multi-TeV squark masses can naturally accommo¬ 
date the stop masses required to achieve a 125 GeV 
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Higgs boson within the MSSM. 

In this paper, we demonstrate that associated 
squark-gaugino production at a y/s = 100 TeV pro¬ 
ton collider provides a probe of > 10 TeV squark 
masses which is complementary to pair production. 
Our main results are summarized in Figures 1-5, 
which show the reach of a ^/s = 100 TeV p-p col¬ 
lider with 3 ab“^ integrated luminosity for squark- 
gaugino associated production in various spectra ^ . 

Squark-gluino production can discover squark 
masses up to 32 TeV for < 4 TeV gluino masses 
in spectra with a large gluino-neutralino LSP mass 
splitting (Fig 1). For spectra with a small gluino- 
neutralino LSP mass splitting, squark masses up to 
37 TeV can similarly be discovered (Fig. 2). No¬ 
tably, our analysis finds that the gluino-neutralino 
DM coannihilation region [18, 19] can be excluded 
for squark masses < 28 TeV. For squark-Wino 
(Bino) LSP production, Wino (Bino) masses up to 
4(1.7) TeV can be discovered for squark masses 
< 7 (5) TeV (Figs. 3-4). We find a similar reach for 
squark-Wino NLSP production (Fig 5), even with¬ 
out utilizing objects resulting from NLSP —>■ LSP 
decay. Our results indicate that squark-gaugino 
production represents a SUSY discovery mode at a 
^/s = 100 TeV p-p collider in a wide variety of mod¬ 
els with heavy first- and second-generation squarks. 

The remainder of this paper is organized as fol¬ 
lows. Section H discusses our general methodology 


^ Note that a recent study in [17] calls for an integrated lu¬ 
minosity of between 10 and 20 ab“^ at a future 100 TeV 
p-p collider. We present here results for 3 ab“^ as a conser¬ 
vative estimate, and so as to be directly comparable with 
the current literature. 




and simulation strategies. Section III presents in 
detail our analysis of squark-gluino associated pro¬ 
duction, while Section IV presents our analysis of 
squark-Wino/Bino associated production. We sum¬ 
marize our results in Section V. 


II. GENERAL METHODOLOGY 

In this section we briefly discuss the general 
methodology of the analyses presented below. Event 
topologies arising from heavy squark - light gaugino 
associated production are characterized by a hard 
leading jet and significant These objects result 
primarily from the squark decay products, as the as¬ 
sociated gaugino is produced at relatively low trans¬ 
verse momentum. The dominant SM background 
for such events is in the tt+ jets and vector boson -|- 
jets channels [ 1 ], which fall off rapidly with increas¬ 
ing leading jet pt, $t-, and [Ht is defined 

as the scalar sum of the jet transverse energies). 

In the following analyses, we consider the reach 
of a y/s — 100 TeV proton-proton collider given 3 
ab“^ integrated luminosity. The minimum produc¬ 
tion cross section yielding > 10 events is roughly 
~ 10“^ fb, corresponding to -I- rrig ~ 35 TeV 
{rriq + ^ 15 TeV) for squark-gluino (squark- 

Wino) associated production. For such masses, good 
background discrimination is achieved with hard 
leading jet pt cuts for squark-gluino production, and 
with hard $x/\/Ht cuts for squark-Wino/Bino pro¬ 
duction. Our strategy is as follows: for each analy¬ 
sis we impose a set of baseline cuts catered to a set 
of spectra. We then scan over leading jet pt and 
cuts (squark-gluino) or $j’/\/Ht cuts (squark- 
Wino/Bino) to maximize significance cr, defined by 


^/l + B + + 72S'2 ■ 

S (B) is the number of signal (background) events 
passing cuts, and 7 (A) parameterize systematic un¬ 
certainties associated with signal (background) nor¬ 
malization. Details of the event generation and col¬ 
lider simulation are given in Appendix A. Like most 
future collider studies, our simulated a values are 
subject to 0 ( 1 ) uncertainties associated with e.g. 
the performance of a detector which is yet to be de¬ 
signed. However, this translates to a comparatively 
mild uncertainty for the predicted reach, due to the 
rapid falling of production cross sections with in¬ 
creasing mass. 


Simplified Models 


In the analyses presented below, we consider the 
following SUSY simplified models: 


Model 

Particle Content 

Fig. 

Squark-Gluino 

N on-compressed 
Compressed 

9, g,Xi= B 

Ml = 100 GeV 
ms — m^o = 15 GeV 

Fig. 1 
Fig. 2 

Squark-Wino LSP 

q, Xi = W 

Fig. 3 

Squark-Bino LSP 

q, Xi^B 

Fig. 4 

Squark-Wino NLSP 

Split 

Non-split 

q, NLSP = W,Xi= B/H 
Ml/ g — 100 GeV 
— m^o = 200 GeV 

Fig. 5 


TABLE 1. Simplified models considered in this paper. 

which encompass a wide array of potential event 
topologies arising from squark-gaugino production. 
We take degenerate first and second generation 
squark masses, and decouple all sparticles not listed 
in Table I. For the squark-gluino non-compressed 
model, our results are not sensitive to the choice 
of Ml = 100 GeV as the LSP is effectively mass¬ 
less for m^o <C rrig. The squark-gluino compressed 
model is motivated by the gluino-neutralino coanni¬ 
hilation region [18, 19]. We choose nug — m^o = 15 
GeV as a fiducial value, though the leading jet px- 
based analysis presented below is robust as long as 
rrig — m^o rrig. For the Wino NLSP models, we 
choose two spectra with differing LSP masses to il¬ 
lustrate the effects of increasing the NLSP-LSP mass 
splitting. In the “non-split” case, we have chosen an 
NLSP-LSP mass splitting of 200 GeV so that the 
NLSP decays to the LSP -I- on-shell SM bosons. 


III. SQUARK-GLUINO ASSOCIATED 
PRODUCTION 

In this section we discuss squark-gluino associated 
production. As this process only involves it can 
be important at a -^/s = 100 TeV p-p collider even 
if rriq + rrig > 35 TeV. If a heavy squark of order 
tens of TeV is produced in association with a gluino 
of mass < 10 TeV, the leading jet from the squark 
decay will be very hard, px ~ mql2. Furthermore 
the neutralino resulting from the decay chain q 
qg ^ 3 qx^ will be very boosted, resulting in large 
$x- These kinematic features result in a striking 
collider signature with very low SM background. 
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We explore the reach in squark-gluino production 
at a ^/s = 100 TeV p-p collider for the two types of 
squark-gluino spectra listed in Table I. For simplic¬ 
ity we assume the LSP is a Bino, and all other neu- 
tralinos/charginos are decoupled. Relaxing this as¬ 
sumption allows squark decays to intermediate neu- 
tralinos/charginos, resulting in additional final state 
objects which can be used for background discrimi¬ 
nation. 



FIG. 1. Experimental reach for squark-gluino associated 
production at a 100 TeV proton collider with 3 ab“^ 
integrated luminosity, for spectra with a ~ 100 GeV LSP 
mass. The solid, long dashed and short dashed lines are 
for and 5, 10, 15% systematic uncertainty for the signal 
respectively. Blue lines indicate 5a discovery reach and 
red lines indicate 95% exclusion limits. We assume 20% 
systematic uncertainty in the background. 

For both non-compressed and compressed spectra, 
we impose the following baseline cuts: 

Ht > 10 TeV, > 20 GeV^^^ 

while for the non-compressed spectra we impose the 
additional cut: 

8 jets with pt > 50 (150) GeV 

The softer cut is optimized for heavier squarks and 
lighter gluinos, while the harder cut is optimized for 
lighter squarks and heavier gluinos. Upon imposing 
these baseline cuts, we then scan over leading jet 
Pt and cuts in order to maximize significance 
cr as defined in (1). We have verified that the opti¬ 
mal cuts render any “background” from gluino pair 
production subdominant to the SM background. 

The results of this analysis are depicted in Figs. 
1 and 2, which show the reach of a ^/s = 100 TeV 



FIG. 2. Experimental reach for squark-gluino associated 
production at a 100 TeV proton collider with 3 ab“^ inte¬ 
grated luminosity for spectra with rrig — m^o = 15 GeV. 
The different lines follow the conventions of Fig. 1. We 
assume 20% systematic uncertainty in the background. 

proton collider with 3 ab“^ of integrated luminos¬ 
ity. The solid, long dashed and short dashed lines 
correspond respectively to systematic uncertainties 
of 5, 10 and 15% for the signal normalization, while 
the background systematic uncertainty is fixed to 
20%. The projected reach is fairly insensitive to 
background systematic uncertainties, as the number 
of background events is quite low due to the hard 
leading jet pr and cuts. 

As is evident from Figs. 1 and 2, a =100 TeV 
collider with 3 ab“^ integrated luminosity can begin 
probing much of the “mini-split” parameter space 
for sufficiently low gluino masses. Final states in 
the compressed spectra yield more $t compared to 
the non-compressed spectra, resulting in the greater 
reach depicted in Figure 2. Notably, with 3 ab“^ 
integrated luminosity the entire neutralino-gluino 
coannihilation region (whose upper endpoint lies at 
Trig « TOy « 8 TeV [19]) can be excluded if the 
squark masses are < 28 TeV. 

It is worthwhile to compare Figs. 1 and 2 to pro¬ 
jected reaches for gluino pair production. Our re¬ 
sults for non-compressed spectra have some overlap 
with [1]^, which considered both pair production and 
associated production in similar spectra with squark 


^ A search optimizing over Ht cuts as opposed to leading jet 
Pt cuts was done in [1], For the spectra in Fig. 1, the Ht 
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masses < 24 TeV. The results of [1] indicate that 
gluino pair production will likely be the discovery 
channel for colored superpartners for the spectra in 
Fig. 1 provided rrig < 14 TeV. On the other hand, 
if the gluino and the LSP are nearly degenerate, 
searches for gluino pair production rapidly lose sen¬ 
sitivity [1]. Thus if the gluino and the LSP are nearly 
degenerate as in the gluino-neutralino coannihilation 
scenario, squark-gluino associated production would 
be a potential discovery channel for colored super¬ 
partners. 

IV. SQUARK-WINO AND SQUARK-BINO 
ASSOCIATED PRODUCTION 

In this section we discuss squark-Wino and 
squark-Bino associated production. These chan¬ 
nels are particularly important if squark-gluino as¬ 
sociated production is inaccessible due to a suffi¬ 
ciently heavy gluino mass^. The event topology is 
qualitatively similar to squark-gluino production, as 
the squark will decay to a boosted jet and boosted 
Wino/Bino while the associated Wino/Bino is pro¬ 
duced at relatively low pT- However as noted in 
Section II, associated squark-Wino/Bino produc¬ 
tion probes significantly lighter squark masses than 
squark-gluino production. Consequently, multi-TeV 
leading jet px and $rp cuts are not as effective 
for background discrimination in squark-Wino/Bino 
production. Instead, we find that hard $x/VHt 
cuts are quite effective at reducing the tt+ jets and 
vector boson -|- jets background without rejecting 
too many signal events. 

In order to determine the projected reach for 
squark-Wino/Bino production at a y/s = 100 TeV 
pp collider with 3 ab“^ integrated luminosity, we 
impose the following baseline cuts: 

PT(ji)>2TeV, |;r>3TeV, A</(j, > 0.5 

where the Acj) cut is imposed only on the two lead¬ 
ing jets. We then scan over $x/VHt cuts for each 
spectrum to maximize a as defined in (1). 

Our focus is on spectra listed in Table I where 
at most one of the gaugino/Higgsino mass param¬ 
eters Ml, M 2 , p are < 1 TeV, such that the gauge 


cut based analysis has a 3-5 TeV weaker reach in rriq -\- rrig 
with respect to squark-gluino associated production. 

^ In the MSSM, a gluino which is hierarchically heavier than 
the squarks requires fine-tuning of the soft masses. This can 
be avoided however in a model with Dirac gluinos [20, 21]. 


eigenstates are approximately aligned with the mass 
eigenstates in the neutralino/chargino sectors. We 
omit the “compressed” region rrig — < 1 TeV, 

as in this region the event topology of associated 
squark-Wino/Bino production is similar to squark 
pair production, only with a substantially smaller 
cross section. Assuming a systematic uncertainty of 
10 % for the signal normalization, the results of the 
above analysis for the various spectra in Table I are 
depicted in Figures 3-5. 



m, (GeV) 

FIG. 3. Experimental reach for squark-Wino LSP associ¬ 
ated production at a 100 TeV proton collider with 3 ab“^ 
integrated luminosity. The solid, long dashed and short 
dashed lines are for 1, 2, 3% systematic uncertainty for 
the background respectively. Blue lines indicate 5a dis¬ 
covery reach and red lines indicate 95% exclusion limits. 
We do not consider the grey shaded region {rriq—m-^ < 1 
TeV) for reasons given in the text. We assume 10% sys¬ 
tematic uncertainty for the signal. 

Figure 3 shows the reach for squark-Wino produc¬ 
tion with a pure Wino LSP; the solid, short-dashed, 
long-dashed lines correspond to background uncer¬ 
tainties of 1%, 2% and 3 %. In Figure 4 we show the 
reach for squark-Bino production with a pure Bino 
LSP. The solid, short-dashed, long-dashed lines cor¬ 
respond to background systematic uncertainties of 
0.5%, 1% and 1.5%. Compared to squark-Wino pro¬ 
duction, the reach for squark-Bino associated pro¬ 
duction is quite sensitive to background uncertain¬ 
ties. This is because the 5a contours for squark-Bino 
production correspond to significantly lower masses 
due to the smaller production cross-section, result¬ 
ing in lower optimal cuts and thus larger 

backgrounds. 

In Figure 5 we show the reach of the 
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FIG. 4. Experimental reach for squark-Bino LSP as¬ 
sociated production at a 100 TeV proton collider with 
3 ab~^ integrated luminosity. The solid, long dashed 
and short dashed lines are for and 0.5, 1, 1.5% system¬ 
atic uncertainty for the background respectively. Blue 
lines indicate 5a discovery reach and red lines indicate 
95% exclusion limits. We do not consider the region 
{niq — rrig < 1 TeV) for reasons given in the text. We 
assume 10% systematic uncertainty in the signal. 


based monojet analysis for squark-Wino production 
with a Wino NLSP, with background uncertainties 
fixed to be 1%. The green lines correspond to 
Mmlsp — Mlsp = 200 GeV, while the red lines 
correspond to Mlsp = 100 GeV. For comparison, 
the blue lines show the reach for squark-Wino pro¬ 
duction when the Wino is the LSP. Away from the 
rriq ~ region the sensitivity is lower for a Wino 
NLSP, as is being traded for W, Z and higgs 
bosons arising from the NLSP —?► LSP decay. Note 
that the analysis considered here does not exploit the 
additional SM bosons present in the Wino NLSP sce¬ 
nario. Thus the reach for the Wino NLSP scenario 
depicted in Figure 5 applies regardless of whether 
the LSP is Bino-like or Higgsino-like. Exploiting the 
additional SM bosons could extend the reach for the 
Wino NLSP scenario, so the result presented here is 
a conservative estimate. 

We close this section by comparing the results of 
Figures 3-5 to studies of pair production at s = 100 
TeV. Given 3 ab“^ integrated luminosity, squark 
pair production can discover squark masses up to 2.5 
TeV [1] (assuming a conservative 20 % background 
systematic uncertainty). In the pure Wino case, 
searches in VBF channels can discover Winos up to 
1.1 TeV [10]. Disappearing tracks can also provide 



mu (GeV) 

FIG. 5. Experimental reach for squark-Wino associated 
production at a 100 TeV proton collider with 3 ab“^ 
integrated luminosity. Solid lines indicate 5a discovery 
reach, and dotted lines indicate 95% exclusion limits. 
Blue curves correspond to a Wino LSP, while the green 
(red) curves correspond to a Wino NLSP with Mnlsp — 
Mlsp = 200 GeV {Mlsp ~ 100 GeV). The results are 
applicable for both Bino- and Higgsino-like LSP. We do 
not consider the grey shaded region (m, — < 1 TeV) 

for reasons given in the text. We assume 1% systematic 
uncertainty in the background and 10% in the signal. 


a collider probe of pure Wino LSP pair production. 
Extrapolating the disappearing tracks background 
from the 8 TeV ATLAS study [22], the projected 
reach is 2-3 TeV for pure Winos [3]. However, the 
data-driven disappearing-track background at 100 
TeV is difficult to estimate, making this projected 
reach less reliable than the reach in the VBF chan¬ 
nel or the reach depicted in Figure 3. Finally, pair 
production of Wino NLSPs has been considered in 
[6, 7]. Assuming no systematic uncertainties, for a 
Higgsino LSP the projected discovery reach is 2.3 
TeV, while for a Bino LSP the reach is 1-3 TeV de¬ 
pending on the NLSP —>■ Z LSP branching ratio. 
Gomparing these reaches to Figures 3-5, we see that 
squark-Wino/Bino associated production can pro¬ 
vide a SUSY discovery mode provided the squark is 
not too much heavier than the Wino/Bino. 


V. SUMMARY 

We have examined in this paper the kinematic 
reach for squark-gaugino associated production at 
a 100 TeV proton proton collider. In models where 
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squark pair production is kinematically inaccessible 
at a 100 TeV collider, squark-gaugino associated pro¬ 
duction may be the discovery mode for SUSY in a 
large portion of parameter space. 

We have considered the various simplified mod¬ 
els listed in Table I. For squark-gluino production 
with C>(TeV) gluinos, the discovery reach for first- 
generation squarks can be up to 37 TeV for com¬ 
pressed spectra (small gluino-LSP mass splitting), 
and up to 32 TeV for non-compressed spectra, sub¬ 
ject to systematic uncertainties. For squark-Wino 
LSP production, we have shown that the discovery 
reach for the Wino is almost 4 TeV for squarks of 
~ 7 TeV, subject to systematic uncertainties. For 
squark-Wino NLSP production we have analysed 
two scenarios: one where the NLSP-LSP mass dif¬ 
ference is 200 GeV, and one where the LSP mass is 
~ 100 GeV. In the first scenario, the Wino discovery 
reach is about 3.5 TeV for squarks of ^ 7 TeV. In 
the second scenario, the Wino reach extends up to 
6 TeV. Our results in the Wino-NLSP scenario are 
insensitive to the nature of the LSP. For < 9 TeV 
squark masses, squark-Wino associated production 
marks a significant increase in the Wino reach com¬ 
pared to pair production channels. We also consider 
squark-Bino associated production, and find that the 
kinematic reach for the Bino is up to 1.7 TeV for 
squarks of mass ^ 5 TeV, subject to systematic un¬ 
certainties. 

The results presented here raise the exciting 
prospect of directly probing a region of parameter 
space that so far has been the exclusive domain of 
indirect searches through low-energy FGNC observ¬ 
ables. The squark-gaugino associated production 
channels studied here, coupled with studies of su¬ 
persymmetry at 100 TeV colliders already under¬ 
taken [1-11], provide a strong physics case for the 


construction of such a collider. 
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Appendix A: Event generation 

Signal events were generated using MADGRAPH5 
[23] , with showering and hadronization implemented 
via PYTHIA6.4 [24]. We do not perform MLM for 
the signal events. We have validated this approxi¬ 
mation by performing MLM with 2 additional jets 
for a number of benchmark spectra. We use the sim¬ 
ulated Snowmass backgrounds [25] , proccessed with 
Delphes3.1.2 [26] supplemented by the Snowmass 
detector card [27] for a ^/s = 100 TeV hadron col¬ 
lider. Production cross sections for squark-gluino 
associated production are computed at NLO using 
PR0SPIN02 [28]. For squark-Wino/Bino production 
we use the LO result computed by MADGRAPH5 . Event 
analysis is performed with MadAnalysis5 [29]. We 
expect our kinematic cuts to effectively remove any 
contamination from QCD backgrounds and pileup 
effects, so we neglect both of these in our analysis. 
Note that for squark-gluino associated production 
in the rriq 3> rrig region, the dijet background may 
not be negligible for non-compressed spectra. For 
these spectra, jet substructure techniques can help 
distinguish signal events from the QGD background 
[30]. 
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